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We report measurements of the decays B"*" —>■ ^'cpK'^ and (l)(f>K° using a sample of 231 

million BB pairs collected with the BABAR detector at the PEP-II asymmetric-energy B Factory 
at the Stanford Linear Accelerator Center. The branching fractions are measured to be B{B'^ 
(j)<f)K+) = (7.5 ± 1.0 (stat) ± 0.7 (syst)) x lO"*^ and B{B° HK°) = (4.lli;^ (stat) ± 0.4 (syst)) 
X 10"^ for a (fxj) invariant mass below 2.85 GeV/c^. 

PACS numbers; 13.20.He, 14.40.Nd 



We report an observation of the decay (f)(l)K^ 
and evidence for 5° (l)(j>K^ along with their corre- 
sponding branching fractions. The decay modes studied 
involve a flavor-changing neutral current b ^ sss transi- 
tion. These charmless transitions can interfere with the 
b CCS process B ^ rjc K, r]c ^ 4>4> and lead to di- 
rect CP violation [H ; the CP asymmetry expected in the 
Standard Model (SM) is zero, so a non-zero CP asym- 
metry would be a sign of new physics. Furthermore, an 
analysis of time-dependent CP violation in cjxpK'^ 
would be sensitive to physics beyond the Standard Model 
and complementary to measurements in the other decays 
that are dominated by the b sss transition. In the SM, 
the partial decay widths for these decays are expected to 
be equal due to the suppression of AI — 1 transitions 
in the electroweak Hamiltonian [2]. Additional interest 
in these final states arises from the possibility of glueball 
production with subsequent decays to cfxp 3] . 

We study the charmless decays B (j)(f>K by work- 
ing below the charm production threshold (m^^ < 2.85 
GeV/c^) to avoid the region dominated by the rjc res- 
onance. Theoretical estimates of these branching frac- 
tions are in the range (1.3 - 4.2) x 10~^ [4, 5] within 
the above kinematic region. The Belle Collaboration has 
previously reported evidence for the decay B~^ (j)(j)K^ 
with a branching fraction of 2.6^J g (stat) ± 0.3 (syst) x 
10~^ for m^^ < 2.85 GeV/c^ [ij; no measurement of the 
branching fraction for i?" (p<j)K^ has previously been 
reported. Throughout this paper, for any given mode, 
the corresponding charge-conjugate mode is also implied. 

The data used in this analysis were collected with the 
BaBAR detector ^ at the PEP-II asymmetric e+e" stor- 
age ring. These data represent an integrated luminosity 
of 209.1 fb^^ collected at a center-of-mass (CM) energy 
\/s — 10.58 GeV, near the peak of the T(4S') resonance, 
plus 21.6 fb~^ collected at a CM energy approximately 
40 MeV below the T(4S'). These are referred to as the on- 
resonance and off-resonance data samples, respectively. 

Charged particles from the e+e" interactions are de- 
tected and their momenta measured by a five-layer, 
double-sided silicon vertex tracker (SVT) and a 40-layer 
drift chamber (DCH) with a helium-based gas mixture, 
placed in a 1.5-T uniform magnetic field produced by a 
superconducting magnet. The charged particles are iden- 
tified using likelihood ratios calculated from the ioniza- 
tion energy loss (dE/dx) measurements in the SVT and 
DCH, and from the observed pattern of Cherenkov light 
in an internally reflecting ring- imaging detector. A K /n 



separation of better than four standard deviations (a) 
is achieved for momenta below 3 GeV/c, smoothly de- 
creasing to 2.5 a at the highest momenta present in the 
S-decay final states. Photons and electrons are identified 
as isolated electromagnetic showers in a CsI(Tl) electro- 
magnetic calorimeter. The detector response is simulated 
with the GEANT4 [1] program. 

The B-meson daughter candidates are reconstructed 
through their decays 4> K^K~ and tt+tt". 
For (f> K~ , we require one charged track to be 

consistent with the kaon hypothesis, the other to be in- 
consistent with the pion hypothesis, and the invariant 
mass to satisfy 1000 < tuk+k- < 1050 MeV/c^. The 
variable mx+x- will be later used in the fit. The K'^ 
candidates are formed from pairs of oppositely charged 
tracks consistent with the pion hypothesis, with a vertex 

probability greater than 0.001 and a reconstructed de- 
cay length greater than 2 mm. We require the invariant 
mass of the two pions to satisfy 486 < m^+„- < 510 
MeV/c2. 

We reconstruct a _B-meson candidate by combining a 
if + or Kg with two 4> candidates. From the kinematics of 
the T{AS) decays, we determine the energy-substituted 
mass mEs = ((Vs/2)^ ~ Pb'^Y^'^ and the energy differ- 
ence AE = Eg — \fs/'2', where p*g and Eg are the recon- 
structed 3-momentum and energy of the B meson calcu- 
lated in the CM frame, respectively, and y/s is the e+e" 
collision energy in the CM. For signal decays, the toes 
distribution peaks near the nominal mass of the B me- 
son and AE peaks at zero. The AE (toes) resolution 
is about 20 MeV (3.0 MeV/c^). We require \AE\ < 0.2 
GeV, Toes > 5.2 GeV/c^, and the invariant mass of the 
pair of (j) meson candidates to be less than 2.85 GeV/c^. 
The average number of reconstructed B candidates per 
event is 1.06 (1.05) for B+ -> #1^+ (B" -> 00^°). In 
events with multiple candidates we arbitrarily select one 
candidate to avoid a potential bias in the shape of the 
variables used in the selection 

Backgrounds arise primarily from random combina- 
tions of tracks in the continuum e^e~ ^ qq {q = u, d, s, c) 
events. Because of the jet-like topology, in contrast to 
the nearly isotropic distribution of final particles from 
the process T(4S') bb, 

the continuum background can be significantly reduced 
by an appropriate choice of variables describing the event 
shape. Discrimination between signal and continuum 
events is obtained using a Fisher discriminant The 
variable combines eleven event-shape variables defined 
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TABLE I; Fitted signal yield, detection efficiency e(%) including tracking, PID efficiency and fit bias correction, daughter 
branching fraction product Yl^i [lUi significance S (a), measured branching fraction B with statistical and systematic uncer- 
tainties for each decay mode. These branching fractions are for m,/,^ < 2.85 GeV/c^. The first uncertainty is statistical, the 
second systematic. 





Mode 


Signal Yield 


e{%) 


UB. (%) 


Sia) 


S(10-'') 


B+ 


^ (j}(j)K+ 


64 ±9 


15.3 


24.2 


12.9 


7.5 ± 1.0 ±0.7 


B° 




1Q+4.1 


12.6 


8.3 


4.2 


4.11^;^ ±0.4 



in the CM frame the polar angles of the B momen- 
tum vector and the B candidate thrust axis with respect 
to the beam axis, and the scalar sum of the momenta of 
charged particles and photons (excluding particles from 
the B candidate) in nine 10° polar- angle intervals coaxial 
with the _B-candidate thrust axis. 

We use Monte Carlo (MC) simulation for an initial es- 
timate of the residual BB background and to identify 
the decays that may survive the candidate selection and 
have characteristics similar to the signal. We find that 
the contributions from the multi-kaon decays, iJ+Z" 
(j)K+K-K+/° and B+/° K+ K' K+ K' , are 
negligible after selecting events with two (j) meson can- 
didates. 

We obtain the signal yields from an unbinned extended 
maximum-likelihood fit. The variables used in the fit are 
AE, mES, the invariant masses of two meson candi- 
dates, and !F. The likelihood function has two categories 
of probability-density functions (PDF), one for signal and 
the other for the continuum background. The likelihood 
function is defined as 



C 



-(En.) 



N 

n 

4=1 



(1) 



where N is the number of candidates, nj is the number 
of events in category j, and Vj('x.i) is the corresponding 
PDF, evaluated with the observables of the «th event. 
Since correlations among the observables are small, we 
take each V as the product of the PDFs for the sepa- 
rate variables. Possible systematic effects arising from 
correlations are discussed later. 

We determine the signal PDF parameters from MC 
simulated data. We generate signal MC assuming that 
the B meson decays isotropically to (l)(t>K, using three- 
body phase space. The signal PDF distributions are 
parametrized using a single Gaussian function for mss, 
a sum of two Gaussian functions with the same mean 
for AE, a sum of an asymmetric Gaussian function with 
a different width below and above its maximum, and a 
single Gaussian for T . The (f) candidate mass distribu- 
tions are parametrized using a relativistic Breit-Wigner 
distribution convolved with a Gaussian resolution func- 
tion. Control samples (e.g., B D{KTnT)7r) are used 
to verify the resolutions obtained from signal MC. The 



signal PDFs are obtained using correctly reconstructed 
B (j>(j)K decays from MC simulated data. 

The background PDFs are determined using ttles and 
AE sideband data (5.20 < toes < 5.26 GeV/c^, 0.1 < 
\AE\ < 0.2 GeV). We use a first-order polynomial for 
AE, an empirical phase-space function ^lOj for toes, and 
an asymmetric Gaussian function for J^. Since the back- 
ground includes both resonant and non-resonant K^K~ 
combinations, the 0-candidate mass distributions are 
parametrized as the sum of the (f> lineshape (as described 
above) and a first-order polynomial. The parameters al- 
lowed to vary in the fit are the signal and background 
yields and all the background PDF parameters except 
the (j) niass and width. The signal yield from a fit per- 
formed on off-resonance data was consistent with zero, 
as expected. 

Before applying the fitting procedure to the data we 
evaluate the possible signal-yield bias from neglecting 
small residual correlations between discriminating vari- 
ables in the signal PDFs. The bias is determined from 
ensembles of mock experiments obtained from samples of 
signal MC events combined with qq background events 
generated from the PDFs. We find a bias of 7% (10%) 
for B+ cj)(l>K+ {B° ~* (j)cj)K^). We correct the signal- 
detection efficiency for this fit bias. 

We compute the branching fractions from the fit- 
ted signal-event yields, detection efficiencies, daughter 
branching fractions, and the number of produced B- 
meson pairs. In Table U we show the fitted signal 
yield, the detection efficiencies, the products of daugh- 
ter branching fractions for each decay mode, the signifi- 
cances S{a), and the measured branching fractions. We 
assume equal decay rates of the T(4S') to B^B~ and 
j^Oj^o rpjj^g statistical uncertainties in the signal yields 
are taken as the change in the central value when the 
quantity — 2 ln£ increases by one unit from its minimum 
value. The significance is taken as the square root of the 
difference between the value of ~21n£ (with systematic 
uncertainties included) for zero signal and its value at 
the minimum. 

In Fig. [1] (a, b), we show the toes projection distribu- 
tions of B~^ (l)(j>K~^ and B° — > (jxpK'^ events with a 
requirement \AE\ < 0.05 GeV. The corresponding AE 
projections for toes > 5.27 GeV/c^ are shown in Fig. [T] 
(c, d). The PDF model represents the data well, and a 
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2.1 2.2 2.3 2.4 2.5 2.6 2.7 2.5 

m,, GeV/c^ 



FIG. 1: The projected rriES distributions of events with 
\AE\ < 0.05 GeV for (a) B+ <^<?!>X+and (b) B° , 
and the projected Ai? distributions of events with rriES > 5.27 
GeV/c^ for (c) B+ (j>(j,K+and (d) B° (l)(j)K'>. Points with 
error bars represent the data, solid lines the total PDF, and 
dashed lines the background PDF. 



FIG. 2: The projected m^^ distributions of events with 
lA^I < 0.05 GeV and tties > 5.27 GeV/c^ for the B+ 
(j)(j)K^ decay mode. The points with error bars represent the 
data in the signal region, and the shaded histogram represents 
the mass distribution of expected background from the AS 
sideband. 



significant signal is seen in i?+ ^ (f)(j)K^ . At the present 
level of statistics, we do not observe any evidence for res- 
onant structure in the (fxpK Dalitz plot. This is consis- 
tent with our use of three-body phase space in the signal 
MC. The invariant mass of two (p mesons from the decay 
5+ (j)(j>K^ is shown in Fig. [51 Both the signal and 
background display smooth behavior with no evidence of 
any structure. We therefore see no evidence to support 
the hypothesis of glueball production. 

The systematic uncertainties are dominated by our 
knowledge of the signal and background PDFs, fit- 
bias correction, signal MC modehng, and possible non- 
resonant background contributions. The PDF-modeling 
error is largely included in the statistical uncertainty 
since most background parameters are free in the fit. The 
uncertainties in the signal PDFs are estimated by vary- 
ing the signal PDF parameters within their errors. We 
estimate the uncertainty to be 3.8% and 4.8% for charged 
and neutral B meson decays, respectively. The system- 
atic uncertainty due to any discrepancy in the signal 
PDFs between the signal MC and the control data sam- 
ples is 1.7% (1.8%) for B+ ^ (l)cj)K+ {B° (l)(l)K°). The 
uncertainty in the fit-bias correction is taken to be a half 
of the correction. To estimate the uncertainty due to the 
non-resonant background, we refit the data by including 
a non-resonant component in the fit. The change in the 
signal yield is taken as a systematic uncertainty; it is 



found to be 5% for the charged B meson decay and 3% 
for the neutral one. The uncertainty due to the use of 
three-body phase space when calculating the signal effi- 
ciency is 3%, as determined by the signal efficiency vari- 
ation across the Dalitz plot. A correction is applied to 
account for known data-MC differences in track-finding 
efficiency. The uncertainty on this correction is 0.8% per 
track. Systematic uncertainty due to the PID require- 
ments are 3.5% and 2.5% for the charged and neutral B 
meson decays, respectively. There is a systematic uncer- 
tainty of 2.1% on the efficiency of iiTg reconstruction. The 
uncertainty on the total number of BB pairs in the data 
sample is 1.1%. Published data [llj provide the uncer- 
tainties in the i?-daughter product branching fractions 
(0.2 - 1.4%). 

In conclusion, in the charged decay mode, we observe a 
signal of 64±9 (stat) events with a significance of 12.9 cr, 
corresponding to a branching fraction of B{B'^ (p(j)K^) 
= (7.5 ± 1.0 (stat) ± 0.7 (syst)) xlO"^ where m^^ < 
2.85 GeV/c^. This result is larger than the previous 
measurement reported by the Belle Collaboration and 
is also larger than theoretical predictions. The decay 
_B+ 4>(f)K^ is not dominated by a narrow glueball state 
with mass below 2.85 GeV/c^. In the neutral mode, we 
observe a signal of lO.Ots^ii (stat) events with a signifi- 
cance of 4.2 (7, corresponding to a branching fraction of 
B{B^ cf)(j)K°) = (4.lti:^ (stat) ± 0.4 (syst)) xlQ-^, 
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where < 2.85 GeV/c^. This is the first evidence 

for the process cpcjiK'^. The decay widths of the 

charged and neutral modes differ by less than 2 a. The 
fact that the observed charmless m^0 spectrum appears 
to extend into the region of the r/c resonance opens the 
possibility of looking for direct CP violation in interfer- 
ence between the two processes. 
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